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employed in the analysis of jets at the LHC and show the a a (Mz)-dependence of the transverse 
EEC and the AEEC in the anticipated range 0.11 < a s (Mz) < 0.13. 
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I. INTRODUCTION 



Hadron jets are powerful quantitative tools to study Quantum Chromo Dynamics (QCD) in high energy physics. 
In e + e~ colliders PETRA, PEP and LEP, and also in the electron-proton collider HERA, jet studies have been 
undertaken extensively. These include the measurements of the inclusive variables, such as thrust, acoplanarity and 
hadron energy flow, as well as the exclusive jet distributions, yielding a consistent and precise value of the QCD 
coupling constant a s {Mz) jH ■ At the hadron colliders Tevatron and the LHC, QCD predictions for jets have been 
compared with the measured transverse momentum (pt) distributions, and also with the multi-jet rates 2 4] assuming 
a jet algorithm JB-0). The theoretical framework for calculating the jet cross sections in hadronic collisions in the 
next-to-leading order (NLO) accuracy has been in place for well over a decade which has been employed in the 

QCD-bascd analysis of the data. 

In comparison to the e + e~ and the ep experiments, event shape variables have so far received less attention in the 
analysis of the data from the hadron colliders, though first results have been lately published on the measurement of the 
transverse thrust and the thrust minor distributions 10( by the CDF collaboration ll| . Studies of the hadronic event 



shapes in pp collisions at the LHC have just started, initiated by the CMS collaboration using the central transverse 
thrust and the central transverse minor variables, where the term central refers to the jets in the central region of 
the detector 12j . This is followed by a similar analysis by the ATLAS collaboration 13[ . The distributions in these 



variables have been compared with a number of Monte Carlo (MC) simulations, with PYTHIA6 jl4j |. PYTHIA8 [15 1 
and HERWIG++ [lij providing a good description of the data. However, a bench-mark in this field, namely a 
quantitative determination of a s (Mz) at the LHC from the analysis of data on event shapes, is still very much a work 
in progress. 

In this paper, we calculate the transverse energy energy correlation and its asymmetry proposed some time ago 17 1 
as a quantitative measure of perturbativc QCD in hadronic collisions. The analogous energy-energy correlation (EEC) 
function measurements - the energy weig hted angular distributions of the produced hadron pairs in e + e~ annihilation 



were proposed by Basham et al. [18[. The EEC and its asymmetry (AEEC) were subsequently calculated in 



0(oq) la, |2flj, and their measurements have impacted significantly on the precision tests of perturbativc QCD and in 
the determination of a s in e + e~ annihilation experiments (for a recent review, see 2l|). Transverse EEC distributions 
in hadronic collisions [It}, on the other hand, are handicapped due to the absence of the NLO perturbative QCD 
corrections. In the leading order in a s (/x), these distributions show marked sensitivities on the renormalization and 
factorization scales fi = fiji and p = pp, respectively, thereby hindering a determination of a s (Mz)- We aim at 
remedying this drawback by presenting a calculation of the transverse EEC function and its asymmetry in 0(ai?(/i)), 
which reduces the scale-dependence to a few per cent. 

The paper is organized as follows. Sec. |H] collects the definitions and some leading-order features of the transverse 
energy-energy correlation. In Sec. IIII| we present the numerical results calculated at next-to-leading order in a s and 
demonstrate that the transverse EEC and its asymmetry are robust against variations of various parameters except 
for a s , for which we present the NLO results in the range 0.11 < a s (mz) < 0.13 at the LHC (-^/i = 7TeV). We 
conclude in the last section. 



II. TRANSVERSE ENERGY-ENERGY CORRELATION AND ITS ASYMMETRY 



We start by recalling the definition of the transverse EEC function [17 1 
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The first sum on the right-hand side in the second of the above equations is over the events A with total transverse 
energy E^; = ^ > E™ m , with the E™ 11 set by the experimental setup. The second sum is over the pairs of 

partons (a, b) whose transverse momenta have relative azimuthal angle cf> to <j> + Atfi. In addition, the fiducial volume 
is restricted by the experimental acceptance in the rapidity variable rj. 

In leading order QCD, the transverse energy spectrum dcr/effix is a convolution of the parton distribution functions 
(PDFs) with the 2 — > 2 hard scattering partonic sub-processes. Away from the end-points, i.e., for 0^0° and 
4> 7^ 180°, in the leading order in a s , the energy- weighted cross section d 2 T, / dE^ dcj> involves the convolution of the 
PDFs with the 2 — > 3 sub-processes, such as gg — > ggg. Thus, schematically, the leading contribution for the transverse 
EEC function is calculated from the following expression: 



1 dY' _ S ai A/«i/p(si)/ aa /pfa)*S aia2 ^ l63&3 
a' # Z aiM f ai/p { Xl )f a2/p (x 2 )*a a ^ b ^ 
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where Y J aia2 ^ >hlh ' 2b3 is the transverse energy-energy weighted partonic cross section, Xi (i = 1,2) are the fractional 
longitudinal momenta carried by the partons, f ai /p{xi) and f a2 /p{x2) are the PDFs, and the * denotes a convolution 
over the appropriate variables. The function defined in Eq. ^ depends not only on <fi, but also on the ratio E T " n / y/s 
and rapidity rj. In general, the numerator and the denominator in Eq. (J5|) have a different dependence on these 
variables, as the PDFs are weighted differently. However, as already observed in 17[, certain normalized distributions 
for the various sub-processes contributing to the 2 — > 3 hard scatterings are similar, and the same combination of 
PDFs enters in the 2 — » 2 and 2 — » 3 cross sections; hence the transverse EEC cross section is to a good approximation 
independent of the PDFs (see, Fig. 1 in |17|). Thus, for a fixed rapidity range |?7| < r/ c and the variable Et/-\/S; one 



has an approximate factorized result, which in the LO in a s reads as 
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In the above equation, n/ is the active quark flavor number at the scale fj, and the hadronization scale A is determined 
by the input a s (mz)- The function F{(f>) and the corresponding transverse EEC asymmetry defined as 
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were worked out in [Tjj in the leading order of a s for the CERN SPS pp collider at yfs = 540 GeV. In particular, it 
was shown that the transverse EEC functions for the gg-, gq- and gg-scattcrings had very similar shapes, and their 
relative contributions were found consistent to a good approximation with the ratio of the corresponding color factors 
1:4/9:16/81 for the gg, gq{= gq) and qq initial states over a large range of <f>. 



III. NEXT-TO-LEADING ORDER RESULTS FOR THE TRANSVERSE EEC AND ITS ASYMMETRY 



We have used the existing program NLOJET++ Q, which has been checked in a number of independent NLO 
jet calculations [22j], to compute the transverse EEC and its asymmetry AEEC in the NLO accuracy for the LHC 
proton-proton center-of-mass energy ^/s = 7 TeV. Schematically, this entails the calculations of the 2 — > 3 partonic 
sub-processes in the NLO accuracy and of the 2—^4 partonic processes in the leading order in a s (/■*)> which contribute 
to the numerator on the r.h.s. of Eq. ([2]). We have restricted the azimuthal angle range by cutting out regions near 
= 0° and <f> = 180°. This would, in particular, remove the self-correlations (a = b) and frees us from calculating the 
0(a 2 ) (or two-loop) virtual corrections to the 2—^2 processes. Thus, with the azimuthal angle cut, the numerator in 
Eq. @ is calculated from the 2—^3 and 2 — > 4 processes to O(a^). The denominator in Eq. @ includes the 2^2 
and 2 — > 3 processes, which arc calculated up to and including the 0(a)?) corrections. 

In the NLO accuracy, one can express the EEC cross section as 
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It is customary to lump the NLO corrections in a so-called if-factor (which, as shown here, is a non-trivial function 
of 4>), defined as 



o 



K EEC (<p) = 1 + ^^G(0) . (7) 
The transverse EEC asymmetry in the NLO accuracy is likewise defined as 

(8) 



a' dej) 

and the corresponding K-factor is defined as 



1 effi— cM 
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The principal result of this paper is the calculation of the NLO functions K EEG {<j)) and K AEEC (<fi) and in demonstrating 
the insensitivity of the EEC and the AEEC functions, calculated to NLO accuracy, to the various intrinsic parametric 
and the underlying event uncertainties. 

We now give details of the computations: In transcribing the NLOJET++ Q program, we have replaced the default 



241 sets. 



structure functions therein by the state of the art PDFs, for which we use the MSTW [23[ and the CT10 
We have also replaced the kr jet algorithm by the anti-fcy jet algorithm for defining the jets, in which the distance 
measures of partons are given by 



d i:j =mm{k ti ,k tj ) — 2 , d. lB = k ti , (10) 



-2 ,-2 



with R being the usual radius parameter. We recall that the NLO corrections we are using [9( have been computed 
in the Catani-Seymour dipole formalism In particular, it involves a certain cutting of the phase space of the 
dipole subtraction terms and the numerical calculations require the generation of a very large number of events (we 
have generated O(10 10 ) events on the DESY-Thcory PC Cluster) to bring the statistical accuracy in the NLO EEC 
distribution to the desired level of below a few per cent. We have assumed the rapidity range I77I < 2.5, have put a cut 
on the transverse energy Et > 25 GeV for each jet and require Eti + Et2 > 500 GeV for the two leading jets. The 
latter cut ensures that the trigger efficiencies for the LHC detectors will be close to 100%. We have set the transverse 
energy of the hardest jet as the default factorization- and renormalization-scale , i.e., /if = = i?™ ax . We then 
vary the scales /if and /i/? independently in the range 0.5i?™ ax < {hf,IJ>r) _ 2 £?™ ax to study numerically the scale 
dependence. 

The effects induced by the un derly ing event, multiparton interactions and hadronization effects have been studied 
by us using the PYTHIA6 MC (l5|. In Fig |TJ we show a comparison of the transverse EEC and its asymmetry 
for R=0.6 and R=0.4 with and without the underlying event effects (UE). In Fig |2j the results of the transverse 
EEC and its asymmetry at the hadron and parton level are presente for R=0.6 and 11=0.4. To better display this, 
we show in Fig. [3] the normalized distribution of the hadronization factor (left) and the underlying events factor 
(right), from which it is easy to see that both the hadronization and the UE effects are small. Typically, the effect of 
hadronization on the transverse EEC is < 5% and from the underlying event < 6% for the jet-size parameter R = 0.6. 
The corresponding numbers are < 5% and < 2% for R = 0.4. The parameter specifying the jet-size in the anti-fc-r 
algorithm is chosen as R — 0.4 in the rest of this paper, as this choice makes the transverse EEC distribution less 
sensitive to the underlying minimum-bias events. Moreover, a smaller value of R induces smaller distortions on the 
EEC distribution for the smaller values of the angle </>. 

An important issue is the effect of the parton showers in the transverse EEC and the AEEC distributions. They 
are crucially important in the 0—^0° and <fi —> 180° angular regions, but their effect is expected to be small in the 
central angular range on which we have concentrated. We have checked this (approximat ely) by comparing the results 



in the LO accuracy with those from the parton shower-based MC generator PYTHIA6 [15j, which is accurate in the 
leading log approximation and also includes some NLO terms. Matching the NLO computations with the parton 
shower simulations in the complete next-to-leading log (NLL) accuracy is the aim of several approaches, such as the 



POWHEG method, pioneered and subsequently developed in [25|, [26j, which would allow to quantitatively compute 



the end-point region in the transverse EEC cross section [271 ] . Likewise, resummed perturbative tec hniq ues have been 



developed in a number of dedicated studies for some event shape variables in hadronic collisions 10|, [28| , which would 
expand the domain of applicability of the perturbative techniques to a wider angular region in </>. 
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FIG. 1. (color online) Differential distribution in coscA of the transverse EEC cross section [(a),(b)] and its asymmetry [(c), (d)] 
obtained with the PYTHIA6 MC program [lj with and without the underlying events at yfs = 7 TeV and the anti-fcr algorithm 
with two assumed values of the jet-size parameter R = 0.6[(a),(c)] and R = 0.4 [(b), (d)]. 



In view of the preceding discussion, we have restricted cos</> in the range [—0.8,0.8] which is sliced into 20 bins for 
the presentation of our numerical results. We first show the dependence of the transverse EEC calculated in the NLO 
accuracy on the PDFs in Fig. @] for the two widely used sets: MSTW J23| and CT10 24 1, using their respective central 
(default) parameters. This figure shows that the PDF-related differences on the transverse EEC are negligible, with 
the largest difference found in some bins amounting to 3%, (but typically they are < 1%). We also remark that the 
intrinsic uncertainties from the MSTW2008 PDFs, taking the first 10 eigenvectors of the PDF sets to evaluate the 
distributions, are found negligibly small in the transverse EEC (at most a few per mill), while in the case of CT10, 
these uncertainties are somewhat larger but still below 1% in the EEC. The insensitivity of the transverse EEC cross 
section to the PDFs provides a direct test of the underlying partonic hard processes. In what follows, we will adopt 
the MSTW [23[ PDF set as it provides a correlated range of a s (Mz) and the structure functions for the current range 
of interest for a s (M z ): 0.11 < a s (M z ) < 0.13. 

We next explore the dependences of the transverse EEC cross section and its asymmetry on the factorization and 
the renormalization scales in the range (hf, Hr) = [0.5, 2] x i5™ ax and display them in Fig. [5jfor the transverse EEC 
and Fig.|B]for the asymmetric transverse EEC. Effects of the variations in the scales hf and /j,r on the transverse EEC 
cross section in the LO are shown in Figs. [3J (a), [Sic) and[S] (c), which arc obtained by setting the scales /Up = fXu, 
fixing fxp = -E™ ax and varying and fixing ji^ ~ and varying ^.p, respectively. The corresponding asymmetry 
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FIG. 2. (color online) Differential distribution in cos<^> of the transverse EEC cross section [(a),(b)] and its asymmetry [(c), (d)] 
obtained with the PYTHIA6 MC program [U at the hadron and parton level at ^fs — 7 TeV for the indicated values of the 
jet-size parameter R. 



of the transverse EEC cross sections are displayed in Figs. [5] (a), |5] (c) and [5] (e). We note that the dominant scale 
dependence in the LO arises from the variation of the renormalization scale fiR. This is understandable as the LO 
matrix element have no ^-compensating contribution. The results obtained in the NLO are shown in FigsJ5](b),[5](d) 
and[5](f) for the transverse EEC and in Figs. [6] (b) , [6] (d) and|6](f) for the asymmetry. One observes significantly less 
dependence on the scales; in particular the marked /Xfl-dependence in the LO is now reduced. Typical scale- variance 
on the transverse EEC distribution in the NLO is found to be 2% - 3%, with the largest effects in some bins reaching 
5%. This scale-insensitivity in the NLO accuracy is crucial to undertake a quantitative determination of a s from the 
collider jet data. 

Having shown that the uncertainties due to underlying events and the PDFs are negligible, and the scale dependence 
is much reduced in the NLO, we present our results for the transverse EEC in the LO and the NLO accuracy in 
Fig. Eta), and the corresponding results for the transverse AEEC in Fig. [7] (b). We also compute these distributions 
from a MC-based model which has the LO matrix elements and multiparton showers encoded. To be specific, we have 
used the PYTHIA8 MC program and have generated the transverse EEC and the AEEC distributions, which 
are also shown in Fig. [3a) and Fig. [T^b), respectively. This comparison provides a practically convenient way to 
correct the PYTHIA8 MC-based theoretical distributions, often used in the analysis of the hadron collider data, due 
to the NLO effects. In Fig. [7] (c), we show the function if EEC (0) defined in Eq. © (denoted as NLO/LO in the 
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FIG. 3. (color online) Normalized distribution of the hadronization factor (left) and underlying events effects (right) obtained 
with the PYTHIA6 MC program for the two indicated values of R. 
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FIG. 4. (color online) Dependence of the transverse EEC cross section (a) and its asymmetry (b) on the PDFs at NLO in a s . 
Red entries correspond to the MSTW PDFs and the black ones are calculated using the CT10 PDF set Q. The errors 
shown reflect the intrinsic parametric uncertainties in each PDF set and the Monte Carlo integration uncertainties. 



figure) and another phenomenological function in which the NLO transverse EEC distribution is normalized to the one 
generated by the PYTHIA8 (l5| MC program (denoted as NLO/PYTHIA). The corresponding function K AEEG ((j)), 
defined in Eq. @, is shown in Fig. [7] (d). Here also we show the corresponding phenomenological function in which 
the transverse EEC obtained in NLO is normalized to the ones generated by the PYTHIA MC. We remark that the 
effects of the NLO corrections are discernible, both compared to the LO and PYTHIA8 15(, and they are significant 
in the large-angle region (i.e., for cos<j> < 0). To summarize the NLO effects in the EEC distribution, they reduce the 
scale-dependence, in particular on /x^, and distort the shape of both the EEC and AEEC distributions, providing a 
non-trivial test of the NLO effects. 

Having detailed the intrinsic uncertainties from a number of dominant sources, we now wish to investigate the 
sensitivity of the transverse EEC and the AEEC on a s (Mz). In relating the strong coupling a s (/i) at a certain 
scale relevant for the collider jets, such as /z = i?™ ax , to the benchmark value a s (Mz), we have used the two-loop 
/3-function and the explicit formula for transcribing a s {^) to a s {Mz) can be seen in Eq. (j4|). Results for the transverse 
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FIG. 5. (color online) Dependence of the transverse EEC on the scales /xf, and (ir in LO (a,c,e) and in the NLO (b,d,f) in a 
for the indicated values of the scales. Figs. (a) and (b) are obtained by setting /if = fJ.R and varying it [if = fJ.R — [0.5, 2] x E™ 
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FIG. 6. (color online) Same as Fig. [5] but for the asymmetric transverse EEC. 



EEC and the AEEC are shown in Fig. [8] (a) and Fig. |8] (b), respectively, for the three indicated values of a s (Mz). 
= O.ll(blue), = 0.12 (red), = 0.13 (black). The scale uncertainties are included only in the curve corresponding to 
a s (Mz) = 0.12, as it is close to the current world average a s (Mz) = 0.1184 [2^ and hence our focus on this value. To 
demonstrate the intrinsic errors in the calculations of the transverse EEC and its asymmetry, we show the percentage 
size of the errors in the lower part of Fig. [8] (a) and Fig. [8] (b), respectively, for a s {Mz) = 0.12. Concentrating 
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FIG. 7. (color online) Comparison of the LO computation (red entries), NLO calculation (blue entries) and the PYTHIA 
(black entries) results for the transverse EEC (a) and its asymmetry (b). Fig. (c) displays the function K EEC (<f) involving the 
ratios NLO /LO (red entries) as defined in Eq. J7} and a phenomenological function obtained by replacing the LO results by the 
PYTHIA MC results (black entries). Fig. (d) shows the corresponding function K AEEC (cj>) for the transverse EEC asymmetry 
defined in Eq. ((9} . The errors shown are obtained by adding in quadrature all the uncertainties except the ones from scale 
variations, as described in text. 



first on the transverse EEC, we see that the bin-by-bin errors are typically +2% and —6% (for |cos0| < 0.6), and 
somewhat larger for | cos</>| > 0.6. A part of this error is of statistical origin in our Monte Carlo based theoretical 
calculations and is reducible, in principle, with the help of a more effective importance sampling algorithm in the 
event generation. However, a part of the error is irreducible, given the current theoretical (NLO) precision. This is 
quantified for the normalized integrated transverse EEC X-section over the cos <f> range shown in the figures above, 
which largely removes the statistical (bin-by-bin) error: 
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The computational error on the transverse AEEC is larger, as shown in Fig.|S](b) for a s {Mz) = 0.12. In particular, 
the errors for the last four bins in the AEEC X-section are large due to the intrinsically small value of this cross-section 
as cos0 — > 0. However, in the region —0.8 < cos</> < —0.4, a clear dependence of the differential transverse AEEC on 
a s (Mz) is discernible. This is also displayed for the normalized integrated transverse AEEC X-section given below 
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FIG. 8. (color online) Transverse EEC cross section (a) and its asymmetry (b) with three values of a s (Mz) = 0.11 (blue), 
= 0.12 (red), and = 0.13 (black). The bottom panel of the figures demonstrate the size of errors (red) and deviations with the 
values of a s (Mz) = 0.11 (blue) and = 0.13 (black) from the results evaluated with a s = 0.12. 



(in units of 10 3 ), in which the last four bins contribute very little: 
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Details of the calculations and numerical results for other values of the parameter R, cuts on prmin, and the 
center-of-mass energies for the LHC and the Tevatron will be published elsewhere. 

IV. SUMMARY 

To summarize, we have presented for the first time NLO results for the transverse EEC and its asymmetry for 
jets at the LHC. These distributions are shown to have all the properties that are required for the precision tests of 
perturbative QCD. In particular, they (i) are almost independent of the structure functions, with typical uncertainties 
at 1%, (ii) show weak scale sensitivity; varying the scale from /i = Et/Z to \x = 2Et, the uncertainties are less than 
5% with the current (NLO) theoretical accuracy, (iii) their dependence on modeling the underlying minimum bias 
events for judicious choice of the parameter R is likewise mild, ranging typically from 2% to 5% as one varies from 
R = 0.4 to R = 0.6, and (iv) preserve sensitivity to a s (Mz); varying a s (Mz) = 0.11 to 0.13, the transverse EEC 
(AEEC) cross section changes approximately by 20% (15%), and thus these distributions will prove to be powerful 
techniques for the quantitative study of event shape variables and in the measurement of a 8 {Mz) in hadron colliders. 
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